Integrated genome-wide analysis of expression quantitative trait loci aids interpretation of genomic association studies.

BACKGROUND: Identification of single nucleotide polymorphisms (SNPs) associated with gene expression levels, known as expression quantitative trait loci (eQTLs), may improve understanding of the functional role of phenotype-associated SNPs in genome-wide association studies (GWAS). The small sample sizes of some previous eQTL studies have limited their statistical power. We conducted an eQTL investigation of microarray-based gene and exon expression levels in whole blood in a cohort of 5257 individuals, exceeding the single cohort size of previous studies by more than a factor of 2. RESULTS: We detected over 19,000 independent lead cis-eQTLs and over 6000 independent lead trans-eQTLs, targeting over 10,000 gene targets (eGenes), with a false discovery rate (FDR) < 5%. Of previously published significant GWAS SNPs, 48% are identified to be significant eQTLs in our study. Some trans-eQTLs point toward novel mechanistic explanations for the association of the SNP with the GWAS-related phenotype. We also identify 59 distinct blocks or clusters of trans-eQTLs, each targeting the expression of sets of six to 229 distinct trans-eGenes. Ten of these sets of target genes are significantly enriched for microRNA targets (FDR < 5%). Many of these clusters are associated in GWAS with multiple phenotypes. CONCLUSIONS: These findings provide insights into the molecular regulatory patterns involved in human physiology and pathophysiology. We illustrate the value of our eQTL database in the context of a recent GWAS meta-analysis of coronary artery disease and provide a list of targeted eGenes for 21 of 58 GWAS loci.

Gene expression markers of age-related inflammation in two human cohorts.

INTRODUCTION: Chronically elevated circulating inflammatory markers are common in older persons but mechanisms are unclear. Many blood transcripts (>800 genes) are associated with interleukin-6 protein levels (IL6) independent of age. We aimed to identify gene transcripts statistically mediating, as drivers or responders, the increasing levels of IL6 protein in blood at older ages. METHODS: Blood derived in-vivo RNA from the Framingham Heart Study (FHS, n=2422, ages 40-92 yrs) and InCHIANTI study (n=694, ages 30-104 yrs), with Affymetrix and Illumina expression arrays respectively (>17,000 genes tested), were tested for statistical mediation of the age-IL6 association using resampling techniques, adjusted for confounders and multiple testing. RESULTS: In FHS, IL6 expression was not associated with IL6 protein levels in blood. 102 genes (0.6% of 17,324 expressed) statistically mediated the age-IL6 association of which 25 replicated in InCHIANTI (including 5 of the 10 largest effect genes). The largest effect gene (SLC4A10, coding for NCBE, a sodium bicarbonate transporter) mediated 19% (adjusted CI 8.9 to 34.1%) and replicated by PCR in InCHIANTI (n=194, 35.6% mediated, p=0.01). Other replicated mediators included PRF1 (perforin, a cytolytic protein in cytotoxic T lymphocytes and NK cells) and IL1B (Interleukin 1 beta): few other cytokines were significant mediators. CONCLUSIONS: This transcriptome-wide study on human blood identified a small distinct set of genes that statistically mediate the age-IL6 association. Findings are robust across two cohorts and different expression technologies. Raised IL6 levels may not derive from circulating white cells in age related inflammation.

A meta-analysis of gene expression signatures of blood pressure and hypertension.

Genome-wide association studies (GWAS) have uncovered numerous genetic variants (SNPs) that are associated with blood pressure (BP). Genetic variants may lead to BP changes by acting on intermediate molecular phenotypes such as coded protein sequence or gene expression, which in turn affect BP variability. Therefore, characterizing genes whose expression is associated with BP may reveal cellular processes involved in BP regulation and uncover how transcripts mediate genetic and environmental effects on BP variability. A meta-analysis of results from six studies of global gene expression profiles of BP and hypertension in whole blood was performed in 7017 individuals who were not receiving antihypertensive drug treatment. We identified 34 genes that were differentially expressed in relation to BP (Bonferroni-corrected p<0.05). Among these genes, FOS and PTGS2 have been previously reported to be involved in BP-related processes; the others are novel. The top BP signature genes in aggregate explain 5%-9% of inter-individual variance in BP. Of note, rs3184504 in SH2B3, which was also reported in GWAS to be associated with BP, was found to be a trans regulator of the expression of 6 of the transcripts we found to be associated with BP (FOS, MYADM, PP1R15A, TAGAP, S100A10, and FGBP2). Gene set enrichment analysis suggested that the BP-related global gene expression changes include genes involved in inflammatory response and apoptosis pathways. Our study provides new insights into molecular mechanisms underlying BP regulation, and suggests novel transcriptomic markers for the treatment and prevention of hypertension.

Whole blood gene expression and interleukin-6 levels.

BACKGROUND: Circulating interleukin-6 levels increase with advancing age and are a risk factor for various diseases and mortality. The characterization of gene expression profiles associated with interleukin-6 levels might suggest important molecular events underlying its regulation. METHODS AND RESULTS: We studied the association of transcriptional profiles with interleukin-6 levels in 2422 participants from the Framingham Heart Study Offspring Cohort using Affymetrix Human Exon 1.0 ST Array. We identified 4139 genes that were significantly associated with interleukin-6 levels (FDR<0.05) after adjusting for age, sex and blood cell components. We then replicated 807 genes in the InCHIANTI study with 694 participants. Many of the top genes are involved in inflammation-related pathways or erythrocyte function, including JAK/Stat signaling pathway and interleukin-10 signaling pathway. CONCLUSION: We identified and replicated 807 genes that were associated with circulating interleukin-6 levels. Future characterization of interleukin-6 regulation networks may facilitate the identification of additional potential targets for treating inflammation-related diseases.

Biomarker discovery in serum from patients with carotid atherosclerosis.

BACKGROUND: Blood-based biomarkers of atherosclerosis have been used to identify patients at high risk for developing stroke. We hypothesized that patients with carotid artery disease would have a distinctive proteomic signature in blood as compared to a healthy control population without carotid artery disease. In order to discover protein biomarkers associated with increased atherosclerotic risk, we used two different strategies to identify biomarkers from patients with clinically defined atherosclerosis who were undergoing endarterectomy for atherosclerotic carotid artery disease. These patients were compared with healthy matched controls. METHODS: Serum was obtained from patients undergoing endarterectomy (EA; n = 38) and compared to a group of age-matched healthy controls (n = 40). Serum was fractionated using anion exchange chromatography and three different surface-enhanced laser desorption/ionization (SELDI) chip surfaces and then evaluated with mass spectrometry (MS) and two-dimensional difference gel electrophoresis (2D-DIGE). RESULTS: A random forest (RF) analysis of the SELDI-MS protein peak data distinguished these two groups with 69.2% sensitivity and 73.2% specificity. Four unique SELDI peaks (4.2, 4.4, 16.7 and 28 kDa, all p< 0.01) showed the greatest influence in the RF model. The EA patients with a history of prior clinical atherosclerotic plaque rupture manifested as either stroke or transient ischemic attack (symptomatic; n = 16) were compared to patients with carotid atherosclerosis but no clinical evidence of plaque rupture (asymptomatic; n = 22). Analysis of the SELDI spectra did not separate these two patient subgroups. A subgroup analysis using 2D-DIGE images obtained from albumin-depleted serum comparing symptomatic (n = 10) to asymptomatic EA patients (n = 10) found 4 proteins that were differentially expressed (p < 0.01) in the symptomatic patients. These proteins were identified as α(1)-antitrypsin, haptoglobin and vitamin D binding protein that were downregulated and α(2)-glycoprotein precursor that was upregulated in the symptomatic EA group. CONCLUSIONS: SELDI-MS data analysis of fractionated serum suggests that a distinct protein signature exists in patients with carotid atherosclerosis compared to age-matched healthy controls. Identification of 4 proteins in a subset of patients with symptomatic and asymptomatic carotid atherosclerosis suggests that these and other protein biomarkers may assist in identifying high-risk patients with carotid atherosclerosis.

Protein expression profiles distinguish between experimental invasive pulmonary aspergillosis and Pseudomonas pneumonia.

We hypothesized that invasive pulmonary aspergillosis (IPA) may generate a distinctive proteomic signature in plasma and bronchoalveolar lavage (BAL). Proteins in plasma and BAL from two neutropenic rabbit models of IPA and Pseudomonas pneumonia were analyzed by SELDI-TOF MS. Hierarchical clustering analysis of plasma time course spectra demonstrated two clusters of peaks that were differentially regulated between IPA and Pseudomonas pneumonia (57 and 34 peaks, respectively, p<0.001). PCA of plasma proteins demonstrated a time-dependent separation of the two infections. A random forest analysis that ranked the top 30 spectral points distinguished between late Aspergillus and Pseudomonas pneumonias with 100% sensitivity and specificity. Based on spectral data analysis, three proteins were identified using SDS-PAGE and LC/MS and quantified using reverse phase arrays. Differences in the temporal sequence of plasma haptoglobin (p<0.001), apolipoprotein A1 (p<0.001) and transthyretin (p<0.038) were observed between IPA and Pseudomonas pneumonia, as was C-reactive protein (p<0.001). In summary, proteomic analysis of plasma and BAL proteins of experimental Aspergillus and Pseudomonas pneumonias demonstrates unique protein profiles with principal components and spectral regions that are shared in early infection and diverge at later stages of infection. Haptoglobin, apolipoprotein A1, transthyretin, and C-reactive protein are differentially expressed in these infections suggesting important contributions to host defense against IPA.

A canine model of septic shock: balancing animal welfare and scientific relevance.

A shock canine pneumonia model that permitted relief of discomfort with the use of objective criteria was developed and validated. After intrabronchial Staphylococcus aureus challenge, mechanical ventilation, antibiotics, fluids, vasopressors, sedatives, and analgesics were titrated based on algorithms for 96 h. Increasing S. aureus (1 to 8 x 10(9) colony-forming units/kg) produced decreasing survival rates (P = 0.04). From 4 to 96 h, changes in arterial-alveolar oxygen gradients, mean pulmonary artery pressure, IL-1, serum sodium levels, mechanical ventilation, and vasopressor support were ordered based on survival time [acute nonsurvivors (< or =24 h until death, n = 8) > or = subacute nonsurvivors (>24 to 96 h until death, n = 8) > or = survivors (> or =96 h until death, n = 22) (all P < 0.05)]. In the first 12 h, increases in lactate and renal abnormalities were greatest in acute nonsurvivors (all P < 0.05). Compared with survivors, subacute nonsurvivors had greater rises in cytokines and liver enzymes and greater falls in platelets, white cell counts, pH, and urine output from 24 to 96 h (all P < 0.05). Importantly, these changes were not attributable to dosages of sedation, which decreased in nonsurvivors [survivors vs. nonsurvivors: 5.0 +/- 1.0 vs. 3.8 +/- 0.7 ml x h(-1) x (fentanyl/midazolam/ medetomidine)(-1); P = 0.02]. In this model, the pain control regimen did not mask changes in metabolic function and lung injury or the need for more hemodynamic and pulmonary support related to increasing severity of sepsis. The integration into this model of both specific and supportive titrated therapies routinely used in septic patients may provide a more realistic setting to evaluate therapies for sepsis.

Proteomic analysis of inflammatory biomarkers in bronchoalveolar lavage.

To assess markers of lung inflammation, we used SELDI-TOF and 2-DE to study changes in bronchoalveolar lavage (BAL) protein in 33 subjects challenged with local bronchial lung endotoxin and saline and in 11 patients with acute respiratory distress syndrome (ARDS). Differences in the SELDI-TOF spectra were assessed by t-test after baseline subtraction, normalization to total ion current and alignment by m/z calibration. The temporal changes in acute inflammatory BAL (6, 24 and 48 h following endotoxin challenge) on hydrophobic binding chip surfaces revealed the differential presence of proteins of 9, 14, 18 and 28 kDa (all p <0.001) in the inflammatory BAL. This differential pattern was also found in the ARDS BAL. Principal component analysis of the entire SELDI-TOF spectra separated normal BAL, experimental and clinical lung inflammation supporting the notion of a distinctive protein pattern associated with acute lung inflammation. An analysis of the hydrophobic fraction of the inflammatory BAL using 2-DE, identified increased levels of apolipoprotein A1, and S100 calcium-binding proteins A8 and A9 in the inflammatory BAL. This pattern was also found in ARDS BAL after immunoblot analysis. These approaches will be useful to improve current methods of monitoring lung inflammation and to identify new therapeutic targets.

TNF-alpha induced c-IAP1/TRAF2 complex translocation to a Ubc6-containing compartment and TRAF2 ubiquitination.

Signaling through tumor necrosis factor receptor 2 (TNF-R2) results in ubiquitination of TRAF2 by the E3 c-IAP1. In this report, we confirm that TRAF2 translocates to a Triton X-100 (TX)-insoluble compartment upon TNF-R2 engagement. Moreover, TRAF2 ubiquitination occurs in this compartment, from which TRAF2 is degraded in a proteasome-dependent manner. Confocal microscopy demonstrated that the TX-insoluble compartment is perinuclear and co-localizes with endoplasmic reticulum (ER) markers. The ER transmembrane Ubc6 bound to c-IAP1 and served as a cognate E2 for c-IAP1's E3 activity in vitro. Furthermore, Ubc6 co-localized with translocated TRAF2/c-IAP1 in the ER-associated compartment in vivo, and a catalytically inactive Ubc6 mutant inhibited TNF-alpha-induced, TNF-R2-dependent TRAF2 degradation. These results indicate that upon TNF-R2 signaling, translocation of TRAF2 and c-IAP1 to an ER-associated, Ubc6-containing perinuclear compartment is required for the ubiquitination of TRAF2 by c-IAP1. Therefore, the ER plays a key role in the TNF-R-mediated signal transduction cascade by acting as a site of assembly for E2/E3/substrate complexes.

Ubiquitin ligase Smurf1 controls osteoblast activity and bone homeostasis by targeting MEKK2 for degradation.

Bone is constantly resorbed and formed throughout life by coordinated actions of osteoclasts and osteoblasts. Here we show that Smurf1, a HECT domain ubiquitin ligase, has a specific physiological role in suppressing the osteogenic activity of osteoblasts. Smurf1-deficient mice are born normal but exhibit an age-dependent increase of bone mass. The cause of this increase can be traced to enhanced activities of osteoblasts, which become sensitized to bone morphogenesis protein (BMP) in the absence of Smurf1. However, loss of Smurf1 does not affect the canonical Smad-mediated intracellular TGFbeta or BMP signaling; instead, it leads to accumulation of phosphorylated MEKK2 and activation of the downstream JNK signaling cascade. We demonstrate that Smurf1 physically interacts with MEKK2 and promotes the ubiquitination and turnover of MEKK2. These results indicate that Smurf1 negatively regulates osteoblast activity and response to BMP through controlling MEKK2 degradation.

Smurf1 facilitates myogenic differentiation and antagonizes the bone morphogenetic protein-2-induced osteoblast conversion by targeting Smad5 for degradation.

Controlled proteolysis mediated by Smad ubiquitination regulatory factors (Smurfs) plays a crucial role in modulating cellular responses to signaling of the transforming growth factor-beta (TGF-beta) superfamily. However, it is not clear what influences the selectivity of Smurfs in the individual signaling pathway, nor is it clear the biological function of Smurfs in vivo. Using a mouse C2C12 myoblast cell differentiation system, which is subject to control by both TGF-beta and bone morphogenetic protein (BMP), here we examine the role of Smurf1 in myogenic differentiation. We show that increased expression of Smurf1 promotes myogenic differentiation of C2C12 cells and blocks the BMP-induced osteogenic conversion but has no effect on the TGF-beta-induced differentiation arrest. Consistent with an inhibitory role in the BMP signaling pathway, the elevated Smurf1 markedly reduces the level of endogenous Smad5, whereas it leaves unaltered that of Smad2, Smad3, and Smad7, which are components of the TGF-beta pathway. Adding back Smad5 from a different source to the Smurf1-overexpressing cells restores the BMP-mediated osteoblast conversion. Finally, by depletion of the endogenous Smurf1 through small interfering RNA-mediated RNA interference, we demonstrate that Smurf1 is required for the myogenic differentiation of C2C12 cells and plays an important regulatory role in the BMP-2-mediated osteoblast conversion.